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We discuss a novel electromagnetic penguin contribution to the transverse helicity amplitudes in B 
decays to two vector mesons, which is enhanced by two powers of raj/A relative to the standard 
penguin amplitudes. This leads to unique polarization signatures in penguin-dominated decay modes 
such as B —> pK* similar to polarization effects in the radiative decay B —> I\* 7, and offers new 
opportunities to probe the magnitude and chirality of flavour-changing neutral current couplings to 
photons. 
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INTRODUCTION 

Decays of B mesons into two charmless mesons provide 
an abundant source of information on flavour- and CP- 
violating phenomena in the weak interactions of quarks. 
In particular, decays to two vector mesons (B —> VV) 
can shed light on the helicity structure of these interac¬ 
tions through polarization studies. While predicted to 
be fundamentally V-A in the Standard Model (SM), a 
deviation from this expectation cannot currently be ex¬ 
cluded. The first observations of B —> VV decays show 
no anomalies in the helicity structure, but point to a 
reduced amount of longitudinal polarization in penguin- 
dominated decays [1]. This has led to theoretical studies 
that reconsider strong interactions effects in B —> VV 
decays [2-4], or invoke new fundamental interactions [5]. 

Any particular B —> VV decay is characterized by the 
three helicity amplitudes Aq (longitudinal), A_, and A+. 
A quark model or naive factorization analysis [6] leads 
to the expectation that for B, i.e. 6-quark, decay the 
helicity amplitudes are in proportions 

A 0 :A^:A+ = l: — :(—) (1) 

mb V m b J 

with A « 0.5 GeV the strong interaction scale and mb « 
5 GeV the bottom quark mass. This expectation has 
been parametrically (not necessarily numerically) con¬ 
firmed [2] in the framework of QCD factorization, which 
provides a theoretical basis for the heavy-quark expan¬ 
sion of B decays to charmless mesons [7]. The hierarchy 
(1) of helicity amplitudes follows from the V-A structure 
of the standard weak interactions. 

In this Letter we point out and discuss an effect which 
has been neglected in all previous studies of B —> VV, 
but which substantially alters the prediction for polar¬ 
ization observables. The effect is connected with electro¬ 
magnetic penguin transitions, and appears only for neu¬ 
tral vector mesons. It leads to the unique feature that 
the transverse electroweak penguin amplitude is domi¬ 
nated by the electromagnetic dipole operator providing 
a signature similar to polarization in radiative decays 


B —> K* 7 [8], but which is easier to access experimen¬ 
tally. 

The effect in question is related to the two diagrams 
shown in Figure 1. When the vector meson Vi is trans¬ 
versely polarized, there exists a large contribution to the 
decay amplitude due to the small virtuality my 2 of the 
intermediate photon propagator. This is in contrast to 
the case of longitudinal polarization, where the photon 
propagator is canceled, and the amplitude is local on the 
scale mb [9]. The large transverse amplitude is best de¬ 
scribed by a short-distance transition b —> D7 (D = d,s), 
followed by the transition of the low-virtuality photon 
( q 2 -C to 2 ) to the neutral vector meson. We shall per¬ 
form a factorization analysis of the amplitude below. 

The calculation of the diagrams in Figure 1 is straight¬ 
forward. The weak interactions are given in terms of the 
standard effective Hamiltonian [10]. We use the conven¬ 
tions of [11], but generalize the electromagnetic dipole 
operators to include both chiralities 

n eS = % E x l D) E cv“ 7 q ? 7 + ..., (2) 

v p—u,c a ——,+ 

Q^ = ~^Da^(l± l5 )F^b, (3) 

where A = V p b V* D . The ellipses denote other opera¬ 
tors (see [11]). In the SM Cf is suppressed by a factor 
mo/mb, hence Qif is usually neglected. The remain¬ 
ing term is then simply denoted by CV 7 Q7 7 . However, in 
generic extensions of the SM, there is no reason to expect 




FIG. 1: Leading contributions to Aag = ^ w (ViV2) defined in 
the text. 
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a suppression of additional contributions to C^ relative 
to C 7 . The coupling of the photon to the quark elec¬ 
tric charge in V 2 implies that the diagrams of Figure 1 
contribute to the electroweak penguin amplitude in the 
general flavour decomposition of hadronic two-body de¬ 
cay amplitudes. Adopting the a* notation of [9] extended 
to allow for the three helicity amplitudes of B —* VV, the 
new contribution to the transverse electroweak penguin 
amplitudes is 


a <ew(TO = T 


2a, 


3i r 


em r =f /? 


m B m b 


m 


( 4 ) 


Vi 


with C* eff taking into account the effect of quark loop 
diagrams (see Figure 1). R :F is a ratio of tensor to (ax¬ 
ial) vector B —> Vi form factors such that R- equals 1 
in the heavy quark limit [12], while R+ is of order mb/ A. 
We note the large enhancement factor ~ 

(mb/ A) 2 , which implies that the first hierarchy in ( 1 ) 
is inverted, rendering the negative-helicity amplitude A_ 
leading over the longitudinal amplitude Aq in the heavy- 
quark limit. Of course, for real values of mb/mv 2 this 
enhancement is compensated by the small electromag¬ 
netic coupling a em = e 2 /( 47 r). For instance, for neutral 
p mesons, we obtain Aa^~ EW (K*p) « 0.02. This should 
be compared to the uncorrected negative-helicity elec¬ 
troweak penguin amplitude 

«3.ew (K*p) =C 7 + C 9 + ° 8 ^ Cl ° + ... « -0.01, (5) 

and the leading QCD penguin amplitude 

a c 4 ~ ( pK*) = C 4 + + • • ■ « -0.055. ( 6 ) 

The Ci are Wilson coefficients for the various penguin 
operators in the effective Hamiltonian [10], and the el¬ 
lipses denote the 1-loop corrections in QCD factorization 
[4], which we have taken into account in the numerical 
estimates. In the SM the corresponding positive-helicity 
amplitudes are suppressed by about an order of magni¬ 
tude relative to the negative-helicity ones as explained 
above. 

There are strong-interaction corrections to the leading- 
order expression (4) from gluon exchange between the 
quark lines in the second diagram of Figure 1, and also 
through hard interactions with the spectator quark (not 
shown in the Figure) in the B meson. Due to factor¬ 
ization as discussed below, these corrections modify only 
the effective b —> D 7 transition at leading order in the 
expansion in A /mb- They have been computed in next- 
to-leading order in the context of factorization of exclu¬ 
sive radiative B decays [13], and can be incorporated 
by substituting C 7 —> C' 7 (first paper of [13], eq. (62)). 
Turning this argument around, the absolute value of 
Aa^~ EW (K*V 2 ) can be obtained from the branching frac¬ 


tion of B —> A "*7 via 


Aa c -^(K*V 2 ) 
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( 7 ) 
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( 0) 2 


with T (0) « 0.28 a tensor form factor. This results 
in \Aal~ EW (K*p)\ = 0.023, close to the leading-order 
estimate from (4). 

We therefore conclude that the new radiative contri¬ 
bution to the negative-helicity electroweak penguin am¬ 
plitude is at least twice as large (and opposite in sign) as 
was previously assumed. For penguin-dominated b —> s 
transitions it is almost half the size of the leading QCD 
penguin amplitude, and should therefore have visible im¬ 
pact on polarization measurements. In case of new inter¬ 
actions generating C^, the corresponding contribution 
to the positive-helicity amplitude (4) should be observed 
against a very small Standard Model background. 


FACTORIZATION ANALYSIS 

Since the existence of an amplitude violating the power 
counting ( 1 ) may appear surprising, we sketch how this 
amplitude emerges and factorizes in soft-collinear effec¬ 
tive theory (SCET) [14]. The notation and method of the 
following discussion is similar to the one in [15]. After 
integrating out the scale m b , SCET formalizes the inter¬ 
action of the static 5-quark held h v with collinear fields 
for the light-like direction n_, in which meson V\ moves, 
and collinear fields for the light-like direction n+ of meson 
V 2 ■ Let x denote the collinear quark held corresponding 
to V 2 , and let V 2 be the meson that does not pick up the 
spectator quark from the B meson. The leading quark 
bilinears that have non-vanishing overlap with (V 21 are 

X»M lT7s)x, Xfi- 7±(1±75)x- (8) 

The subscript _L denotes projection of a Lorentz vector 
on the plane transverse to the two light-cone vectors nip. 
Both operators scale as A 4 according to the SCET scaling 
rules; the hrst overlaps only with the longitudinal polar¬ 
ization state of V 2 , the second only with a transverse 
vector meson. However, the second operator is not gen¬ 
erated by the V-A interactions of the SM (at least at the 
tree and 1-loop level). This implies the power suppression 
of Azp relative to A 9 in ( 1 ), since the leading contribution 
to transverse polarization now involves an operator with 
an additional derivative D± ~ A 2 ~ A /mb- 

This reasoning ignores electromagnetic effects. Includ¬ 
ing QED in SCET, there is a collinear photon field with 
unsuppressed interactions with collinear quarks (of the 
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same direction). Only the transverse photon field is truly 
a degree of freedom of the theory, since the other two 
components are either gauge-artefacts, or can be elim¬ 
inated by the field equations. Hence there is an addi¬ 
tional operator eA% ± = W^iD t ^ ± W 1 (where W 7 is an 
electromagnetic Wilson line formally required to make 
the operator gauge-invariant), which overlaps only with 
a transversely polarized vector meson. To first order in 
the electromagnetic coupling the matrix element can be 
computed exactly yielding 

(V 2 \[W}iD!; ± W 7 m\0) = ~a V2 ^e ? (9) 

with €j_ a transverse polarization vector, fy 2 the de- 
cay constant, and ay 2 a constant that depends on the 
quark-flavour composition of V 2 , a p = 3/2, a w = 1/2, 
= —1/2. (The convention for the covariant derivative 
corresponding to (3) is iDt/ = id^ ± e q A1/ with e q the 
quark electric charge.) The crucial point is that the op¬ 
erator W^iD^ l± W 7 scales with A 2 , hence this contribution 
to Az p is a factor m b /A larger than even the longitudi¬ 
nal amplitude Aq. Thus, we find the tree-level matching 
equation (see also [16]) 


Q% - [^W!T75 )h v ](0)[W^Wi; ± W 7 ](0), 

( 10 ) 

valid as an equation for the (V 1 V 2 I • • • | B) matrix element. 
In SCET only soft fields can couple to the two brackets 
representing collinear field products in the two different 
directions. But since the photon operator in the second 
bracket is a colour-singlet, the soft fields decouple, and 
the matrix element of the right-hand side of (10) falls 
apart into (9) and (ki|£kP 7 ^(IT 75 ) h v \B), which is pro¬ 
portional to the SCET form factor £j_ [12] at maximal re¬ 
coil. Eq. (10) has to be amended by radiative corrections 
as well as a second operator structure with an additional 
transverse derivative in the first bracket. This is very 
similar to heavy-to-light form factors [15], in fact, these 
corrections simply restore the QCD tensor form factor. 
Combining (9,10), we therefore find 


(ViV^C'?), Q? 7 \B) = im V2 m B 2T 7 1 (0)f V2 a V2 
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(ii) 


which on accounting for the normalization of 0 k^’ew [4, 9] 
reproduces (4). The previous equation should be under¬ 
stood such that the matrix element of Q 7l (Q 7j ) takes 
the value given only when both V\ and V 2 have nega¬ 
tive (positive) helicity, but is zero otherwise. In general, 
the four-quark operators from the effective weak Hamil¬ 
tonian also contribute to the matching coefficient of the 
SCET operator on the right-hand side of (10), and in¬ 
cluding further spectator-scattering effects replaces C 7 
by C' 7 as discussed above. 


THE B -* pK * SYSTEM 


We now focus on the eight B —> pK* decay modes, where 
the electroweak penguin amplitude is largest relative to 
the leading QCD penguin amplitude ( a p = 3/2). Assum¬ 
ing isospin symmetry, the pK* system is described by 
six complex strong interaction parameters for each helic¬ 
ity h = 0,—,+. Neglecting the colour-suppressed elec¬ 
troweak penguin amplitude and the doubly CKM sup¬ 
pressed QCD penguin amplitude is a good approxima¬ 
tion for elucidating the effect of the new (colour-allowed) 
electroweak penguin contribution, hence we write 


A h (p~K*°) 

V2A h (p°K*~) 

A h (p + K*~) 

-V2A h (p°K*°) 


Ph 

[Ph + Ph W \ + e” i7 [T h + C h ] 

Ph + e-^Th 

[P h -P^ w ]+e-^[-C h ], (12) 


and define Xh = Xh/Ph , where Ph is the QCD penguin 
amplitude. The tree amplitudes Th, Ch are suppressed 
by the CKM factor e KM = IKb^l/IWTJ ~ 0.025. 
Assuming 7 = 70° is known, one can obtain Ph from 
an angular analysis of the p~K*° final state, th from 
p±K* z F , and p// n and Ch from the remaining four de¬ 
cay modes. In principle, this allows for a determination 
of P// w , which can be compared to the theoretical re¬ 
sult. In practice, a complete amplitude analysis will be 
experimentally difficult. 

The sensitivity to the electroweak penguin amplitude is 
made apparent in CP-averaged helicity-decay rate ratios 
such as 


2f h (p°K*°) 
f h (p-K*°) 




-p E h W \ 2 + A, 


Ui> 


(13) 


where depends on Ch (and mildly on p// w ), and van¬ 
ishes for Ch —> 0. To estimate £_, we assume that the 
positive-helicity amplitudes are negligible as predicted in 
the SM, and use the observed p~K*° branching fraction 
and longitudinal polarization fraction to determine 
the magnitude of Po and P_. We shall also assume that 
the phase of p// w is not more than 30° away from 0 or 
7 r. Writing pj/ w = [P// w /Th] x th, this amounts to the 
assumption that no large CP asymmetries will be found 
in B —> p ± K* z t 1 . For all other quantities we perform a 
calculation in the QCD factorization framework. In this 
procedure there is a considerable uncertainty in P_ due 
to the discrepant experimental results on / l{p + K *°) [1], 
which may result in an over-estimate of P_ and hence an 
under-estimate of p I / w . It is therefore not excluded that 
the electromagnetic penguin effect is more pronounced 
than in the following theoretical estimates. Keeping this 
in mind, we find R e(p I / w ) = —0.23 ± 0.08 [+0.14lg;Qg] 
and A_ = —0.0 ± 0.2, yielding 


S- = 1.5 ±0.2 [0.7 ±0.1]. 


(14) 
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Here (and below) the numbers in brackets refer to the cal¬ 
culation without the new electromagnetic penguin con¬ 
tribution. Despite the current large theoretical uncer¬ 
tainties, which could be removed with more experimental 
data, eq. (14) clearly shows the impact of this contribu¬ 
tion on polarization observables. The effect is even more 
significant for the ratio of the two final states with neu¬ 
tral p mesons, as S-/S'_ [(15) below] changes by a factor 
of about 4 whether or not the electromagnetic penguin 
contribution is included, but for this ratio the tree con¬ 
tamination is also larger. Data is currently not available 
to test (14), but we may instead consider 


_ 2 r h (p°K*~) 
h f h (p-K*°) 


i +Ph +K- 


(15) 


Following the same strategy as above, we obtain A'_ = 
-0.1 ± 0.0, and S’_ = 0.5 ± 0.1 [1.2 ± 0.1], In the ab¬ 
sence of direct CP asymmetries S' h is directly related to 
the corresponding ratio of polarization fractions f' h = 
fh(p°K*~)/fh(p~K*°). Including a theoretical estimate 
of the CP asymmetries, we obtain 


fo = 1-3 ±0.1 [1.1 ±0.1], (16) 

f = 1 ~^ L [ pR * n ] =0.4±0.1 [0.8±0.1]. (17) 
1 — f L (p~ K*°) 1 1 K ’ 

This can be compared to the experimental values /olexp = 
l-45±8:f|, /ilexp = 0.121°;" [1]. 

Finally we comment on the possibility of detecting 
the presence of new flavour-changing neutral currents in 
the form of an electromagnetic penguin operator with 
opposite chirality, Q^. For this analysis, one must 
isolate experimentally the positive-helicity amplitudes. 
Theoretically, all positive-helicity amplitudes are sup¬ 
pressed, except for the electromagnetic penguin contri¬ 
bution A pEW to the electroweak penguin amplitude. 
In the naive factorization approximation X + = rX_, 
where r is a A/mf,-suppressed form factor ratio, while 
A Pf w « C^/C^AP™ is suppressed only by the 
ratio of Wilson coefficients (see (11)). A conservative 
analysis of the b —> sy branching fraction constrains 
Cf /C^ < 0.5, hence it is possible that the suppres¬ 

sion is weak. This would lead to P^ w P + , in which 
case the positive-helicity decay rates of the p°K* final 
states are much larger than the p ± K* ones. A complete 
angular analysis of the pK* system should allow a deter¬ 
mination of p+ w even when it is not dominant, possibly 
allowing a limit on )C^ of order r « 0.1. 

In conclusion, we discussed an electromagnetic pen¬ 
guin contribution to non-leptonic B decays that has pre¬ 
viously been overlooked. It is the largest contribution 
to the negative-lielicity electroweak penguin amplitude, 
and substantially modifies the theoretical expectations 
for polarization observables in b —> s penguin-dominated 
decays, in particular to the p°K* final states. These ob¬ 
servables may therefore be of considerable interest to the 


search for electromagnetic flavour-changing neutral cur¬ 
rents with chirality equal or opposite to the SM. 
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